MM.shuffle <- MM[sample(1:nrow(MM)),c("Caste","Mass")]
MM.shuffle$Energy <- MM$Energy
fm4 <- lm(Energy ~ Mass + Caste, data = MM.shuffle)
MM.shuffle <- MM.shuffle %>% mutate(pred4 = predict(fm4))
ggplot(MM.shuffle, aes(x = Mass, y = Energy, color = Caste)) +
geom_point(size = 4) +
geom_line(aes(x = Mass, y = pred4, color = Caste), lwd = 2) +
theme(legend.justification = c(0, 1), legend.position = c(0.05, 1)) +
labs(x = "ln Body Mass (g)", y = "ln Daily Energy Expenditure (kJ)") +
scale_color_viridis_d()